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Abstract 

The adjustment of two different selfocs is considered using both exact formulas for the mode- 
connection coefficients expressed in terms of Hermite polynomials of several variables and a qual- 
itative approach based on the Prank-Condon principle. Several examples of the refractive-index 
dependence are studied and illustrative plots for these examples are presented. The connec- 
tion with the tomographic approach to quantum states of a two-dimensional oscillator and the 
Frank-Condon factors is established. 

Keywords: Prank-Condon principle, entanglement, optical waveguide, molecular spectra, Hermite 
polynomials. 

1. Introduction 

It is known that radiation beams propagating in optical waveguides can be described by Schrodinger- 
like equations [1-6]; this has been studied in [7-20]. In this equation, the longitudinal coordinate 
z plays the role of time and the radiation wavelength plays the role of the Planck's constant. The 
refractive-index profile is an analog of the potential energy in the quantum Schrodinger equation. 
The electromagnetic radiation modes in the optical waveguide are analogs of the wave functions of 
a quantum system. This quantumlike picture of the field propagating in the waveguide was widely 
used to study the energy distribution among the modes of successive two (or several) waveguides with 
different refractive index profiles, if the distribution in the first waveguide is known. This problem 



is equivalent to the quantum-mechanical problem of the vibronic structure of electronic lines in the 
absorption and emission spectra of polyatomic molecules. In the case of the harmonic potential, the 
structure is determined in terms of the Frank-Condon factor expressed through Hermite polynomials 
of several variables [21-23] (see, also [24-30]). Recently the technique of femtosecond pulses was 
developed and applied in different domains of ultrapower laser physics [3ll [32] where the problem of 
waveguiding such kinds of radiation also arises [331 El] • 

Our aim here is to consider analogs of the quantum problems of the molecular spectra discussed 
in [27-30] to transfer the results of these investigations into the domain of femtosecond-pulse laser 
physics. The other aim is to connect the classical problem of optical waveguides with the quantum 
problem of entanglement |35l[36]- The point is that for states with two degrees of freedom the notion of 
entanglement corresponds to the degree of correlation of observables related to these different degrees 
of freedom. In the waveguide picture, an analog of the entanglement is related to the structure of 
the modes propagating in the waveguides. When the initial separable two-mode field in the first 
waveguide is propagating into the second waveguide with a different profile of the refractive index, 
the entanglement appears in the arising modes in the second waveguide. The entanglement can be 
considered and related to the energy distribution among the field modes in the second waveguide. 
Here we discuss this relation and consider the possibility to find connections with different criteria of 
the entanglement [37-41]. 

The paper is organized as follows. 

In Sec. 2, we review the results related to the modes in planar waveguides with quadratic refrac- 
tive index. In Sec. 3 the two-mode problem in selfoc-like waveguides is considered, and an analog of 
the two-mode squeezed-light wave function is used for studying the Frank-Condon factor. The en- 
tanglement analogs of the waveguide modes are discussed in Sec. 4. The conclusion and perspectives 
are presented in Sec. 5. 

2. Planar Waveguides 

In this section, we consider the light-beam propagation from one planar waveguide with a quadratic 
refractive index to another waveguide with the propagation direction along the z coordinate. The 
refractive index, being constant along the y and z axes, is described by a symmetric potential-energy 
curve of a one-electron atom as follows: 



The second waveguide has a different dependence of the refractive index, and its axis is shifted from 
the axis of the first waveguide by d 



which is equivalent to shifting and stretching or shrinking of the potential-energy curve. Here we take 
m = m' = 1. 
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Fig. 1. The plot of P^' for J juj = 3 and ujd^ = 9. Fig. 2. The plot of P^' for J juj = 3 and ujd^ = 16. 

The most probable energy level to which the electron jumps can be found in view of the Frank- 
Condon principle. Taking into consideration that the energy levels in a parabolic potential are 
distributed as = + 1/2), the final state can be easily found. 

We consider a one-electron atom. The nuclear Hamiltonian in the Born-Oppenheimer approxi- 
mation reads ^ 

H = -hiu{a, a^}, (3) 

where u is the frequency corresponding to the coordinate q, a and are the creation and annihilation 
operators, respectively. The wave function of electron in the atom with Hamiltonian ([s]) is described 
by 

* = (."■^2"„«(.))-"%xp (-^) H„ (^) . (4) 



where l{uj) = yh/uj, n is a vibrational quantum number, and Hn{^) is the nth Hermite polynomial. 
Shifting the center of the potential from x = to x = d and changing uj ^ uj' and n ^ n' provides 
the wave function 

-1/2 / ~ dY\ TT f X - d^ 



^ {x,n',u',d) = {.y^2-'n'\l{u')y ' exp ( -^^^ ) ( ^TV ) " 
The overlap integral 



2P{uj') j " \ l{uj') 

{n\n') = J ip{x,n,u;,0)'tp{x,n' ,uj' , d)dx (6) 
describes the transition probability amplitude for 

ip{x, n, u, 0) ip{x, n' , u', d). (7) 
The integral ([6]) can be expressed through the Hermite polynomial of two variables 

2l{uj)l{uj') ( X 
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The arguments of the Hermite polynomial read 
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The transition probability from the state \n) to the state \n') is equal to 

Pn = i{n\n')y. (10) 

To illustrate the probability distribution, we present in Fig. 1 a plot of the function P^' for the initial 
level n = 0, potential stretching uj' /uj = 3 and shift ud"^ = 9. The maximum probability is observed 
for the thirteenth level. This result is in agreement with the number calculated directly from the 
Frank-Condon principle. The plot for the transition from n = 3 with u'/u = 3 and ucP = 16 is 
shown in Fig. 2. In this case, the most probable final state is |5). 



3. Elliptic Waveguides 



Now we extend our consideration to the case of elliptical waveguides with quadratic refractive index, 
which corresponds to a one-electron atom with the 2D parabolic potential. 



n{x,y) = k^x'^ + kyU^ ^ U{x,y) 



In complete analogy with the previous section, the shifted and deformed potential reads 



The unshifted wave function of the electron in potential (11) is 

1 / x'^ I 
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In the same way, we can find the probability distribution in the 2D case for the transition 

i) (x, y, n^, Uy, u^, Uy, 0, 0) ^ z/^ (x, y, <, n'y, u'^, u'y, 4, dy^ 
by calculating the overlap integral 

{n\n') = / ip{x,y,n^,ny,u}r,,ujy,0,0)'ilj (^x,y,n'^,n'y,uj'^,ujy,d^,dy'^ dxdy. 
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Figure 3 presents the distribution for the transition from the ground state for u'^/ujx = 2, oj'y/uy 
3, uJxd'^^ = 9, and ujydy = 16. The maximum of this function is observed for the final state |23, 8). 



Fig. 3. Plot of Pq o' ^ for cj^/w^. = 2, to' Juy = 3, Fig. 4. Plot of l' " for cj^Wa; = 2, wi/wj^ = 3 



coxd^ = 9, and Wyd: 



16. 



iOxd'^ = 16, and coydy = 16. 



In Fig. 4, a 3D plot of the transition probability for the initial state |2, 1) and |^ = 2, |^ = 3, 

and uj^dl = ujydy = 16 is presented. The state |8,4) is the most probable final state. 

The pictures are obviously very similar to the planar case. The transition probability distribution 
for the initial state |0, 0) has the Gaussian form. For the transitions from the excited state, one also 
can observe a multi-maxima surface. 



Since there is no xy component in potential (12), the wave created as a result of transition does 



not contain the entanglement term. To obtain a state analogous to the entangled state, one should 
consider the potential with such a term, i.e.. 



U{x,y) 



After the transition, the wave function has the form 



exp 



+ 



+ '^xy 



Hn^ {ax + by) Hn^ {cx + dy) . 
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(n2^.+nynjny\l{ujx)l{ujy)f^ 
As a result, the probability distribution is changed slightly in this case, so we present just plots for 



potential (12) with 7 = 0, keeping in mind that the real plots are very close to these ones. 



4. Conclusions 



To conclude, we point out the main results of the paper. 
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We reviewed the approach for describing the hght propagation in optical waveguides using the 
Fock-Leontovich approximation. In this approximation, the equation describing the hght intensity 
and phase is identical to the Schrodinger equation. In view of this, following [21] we applied an 
analog of the Frank-Condon principle to express the energy distribution among the modes in the 
optical waveguide by the Frank-Condon factors. For harmonic potential, these factors are given 
in the explicit form in terms of special functions associated with Hermite polynomials. We made 
an analysis of some concrete examples and investigated the energy-distribution function for the light 
modes in the optical waveguide. The distribution functions were shown to have maxima corresponding 
to the Frank-Condon principle. 
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